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We extend our recent work on the effects of a time- varying fine-structure constant a in the cosmic 
microwave background, by providing a thorough analysis of the degeneracies between a and the other 
cosmological parameters, and discussing ways to break these with both existing and/or forthcoming 
data. In particular, we present the state-of-the-art CMB constraints on a, through a combined 
analysis of the BOOMERanG, MAXIMA and DASI datasets. We also present a novel discussion of 
the constraints on a coming from large-scale structure observations, focusing in particular on the 
power spectrum from the 2dF survey. Our results are consistent with no variation in a from the 
epoch of recombination to the present day, and restrict any such variation to be less than about 
4%. We show that the forthcoming MAP and Planck experiments will be able to break most of 
the currently existing degeneracies between a and other parameters, and measure a to better than 
percent accuracy. 



I. INTRODUCTION 



The search for observational evidence for time or space 
variations of the 'fundamental' constants that can be 
measured in our four-dimensional world is an extremely 
exciting area of current research, with several indepen- 
dent claims of detections in different contexts emerging 
in the past year or so, together with other improved con- 
straints D, g, §, g §, g, We will review these in Sect. 
II. 

Most of the current efforts have been concentrating 
on the fine-structure constant, a, both due to its obvi- 
ously fundamental role and due to the availability of a se- 
ries of independent methods of measurement. Notewor- 
thy among these is the Cosmic Microwave Background 
(CMB) §, I The latest available CMB results § 
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yield a one-sigma indication of a smaller a in the past, 
but are consistent with no variation at the two-sigma 
level. However, these results are somewhat weakened by 
the existence of various important degeneracies in the 
data, and furthermore not everybody agrees on which 
and how strong these degeneracies are [[)[ [to[ O . 



Here we aim to clarify this issue by analyzing these pos- 
sible degeneracies in some detail, mainly by means of a 
Fisher Matrix Analysis (FMA), see Sect. [TV. We will em- 
phasize that there are crucial differences between 'theo- 
retical' degeneracies (due to simple physical mechanisms) 
and 'experimental' degeneracies (due to the fact that each 
CMB experiment only probes a limited range of scales, 
and that the experimental errors are scale-dependent). 
We will also show how such degeneracies can be elimi- 
nated either by using complementary data sets (such as 
large-scale structure constraints, see Sect. HI) or by ac- 
quiring better data (such as that to be obtained by MAP 
and Planck). We present our conclusions in Sect. |v|. 



In a companion paper fl2|| we will discuss a further 
way in which these degeneracies can be broken, namely 
by including information from CMB polarization. 
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II. THE PRESENT OBSERVATIONAL STATUS 

The recent explosion of interest in the study of varying 
constants is mostly due to the results of Webb and col- 
laborators EL pi 0, II of a 4(j detection of a fine-structure 
constant that was smaller in the past, 

— = (-0.72 ±0.18) x 1CT 5 , z~ 0.5 -3.5; (1) 
a 

indeed, more recent work []l3| provides an even stronger 
detection. These results are obtained through compar- 
isons of various transitions (involving various different 
atoms) in the laboratory and in quasar absorption sys- 
tems, using the fact that the size of the relativistic cor- 
rections goes as (aZ) 2 . A number of tests for possible 
systematic effects have been carried out, all of which have 
been found either not to affect the results or to make the 
detection even stronger if corrected for. 

A somewhat analogous (though simpler) technique 
uses molecular hydrogen transitions in damped Lyman- 
a systems to measure the ratio of the proton and elec- 
tron masses, fi = m p /m e (using the fact that electron 
vibro-rotational lines depend on the reduced mass of the 
molecule, and this dependence is different for different 
transitions). The latest results j?J using two systems at 
redshifts z ~ 2.3 and z ~ 3.0 are 



— = (5.7 ±3.8) x 10~ 5 
A* 



— = (12.5 ±4.5) x 10" 
A* 



(2) 



(3) 



depending on which of the (two) available tables of 'stan- 
dard' laboratory wavelengths is used. This implies a 1.5(7 
detection in the more conservative case, though it also 
casts some doubts on the accuracy of the laboratory re- 
sults, and on the influence of systematic effects ingeneral. 

We should also mention a recent re-analysis M of the 
well-known Oklo bound fll4| . Using new Samarium sam- 
ples collected deeper underground (aiming to minimize 
contamination), these authors again provide two possi- 
ble results for both a and the analogous coupling for the 
strong nuclear force, a s , 



a a s 
a a, 



(0.4 ±0.5) x lO^V -1 



or 



— = -(4.4 ± 0.4) x lO^V -1 ■ 

a s 



(4) 



(5) 



Note that these are given as rates of variation, and effec- 
tively probe timescales corresponding to a cosmological 
redshift of about z ~ 0.1. Unlike the case above, these 
two values correspond to two possible physical branches 
of the solution. See ^ for a discussion of why this 
method yields two solutions (and also note that these 



results have opposite signs relative to previously pub- 
lished ones |f[|). While the first of these branches pro- 
vides a null result, (^J) is a strong detection of an a that 
was larger at z ~ 0.1, that is a relative variation that 
is opposite to Webb's result ([!]). Even though there are 
some hints (coming from the analysis of other Gadolin- 
ium samples) that the first branch is preferred, this is 
by no means settled and further analysis is required to 
verify it. 

Still we can speculate about the possibility that the 
second branch turns out to be the correct one. Indeed 
this would definitely be the most exciting possibility. 
While in itself this wouldn't contradict Webb's results 
(since Oklo probes much smaller redshift and the sug- 
gested magnitude of the variation is smaller than that 
suggested by the quasar data), it would have striking ef- 
fects on the theoretical modelling of such variations. In 
fact, proof that a was once larger than today's value 
would sound the death knell for any theory which mod- 
els the varying a through a scalar field whose behaviour 
is akin to that of a dilaton. Examples include Beken- 



stein's theory [|16| or simple variations thereof 17, [L8| 



Indeed, one can quite easily see 19, 20 that in any such 
model having sensible cosmological parameters and obey- 
ing other standard constraints a must be a monotonically 
increasing function of time. Since these dilatonic-type 
models are arguably the simplest and best-motivated 
models for varying alpha from a particle physics point 
of view, any evidence against them would be extremely 
exciting, since it would point towards the presence of 
significantly different, yet undiscovered physical mecha- 
nisms. 

Finally, we also mention that there have been recent 
proposals Q of more accurate laboratory tests of the 
time independence of a and the ratio of the proton and 
electron masses fi using monolithic resonators, which 
could improve current bounds by an order of magnitude 
or more. 

However, given that there are both theoretical and ex- 
perimental reasons to expect that any recent variations 
will be small, it is important to develop tools allowing 
us to measure a in the early universe, as variations with 
respect to the present value could be much larger then. 

In what follows we focus on the analysis of CMB data 
allowing for possible variations of the fine-structure con- 
stant. In our previous work we have carried out a 
joint analysis using the most recent CMB (BOOMERanG 
and DASI) and big-bang nucleosynthesis (BBN) data, 
finding evidence at the one sigma level for a smaller al- 
pha in the past (at the level of 10 -2 or 10~ 3 ), though at 
the two sigma level the results were consistent with no 
variation. However, as can be seen by comparing with 
earlier work |)|, [nj (and has also been discussed explicitly 
in these papers), these results are quite strongly depen- 
dent on both the observational datasets and the priors 
one uses. 

Regarding this latter issue, we point out that a recent 
p2| improved analysis of standard BBN (focusing mostly 
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on nuclear physics aspects) suggests that 7 Li could lead 
to more stringent constraints on the baryonic density of 
the universe (fib) than deuterium. The point made by 
the authors is that 7 Li is effectively a better baryometer 
than D, because of difficulties in obtaining (extrapolated) 
primordial abundances of the latter. They then obtain 
values for Wb = f2&/i 2 that are considerably lower than 
the standard ones. These results are also corroborated 
by Using these results as a prior would transform 

our previous result ^ into a detection of a varying a at 
more than two sigma. 

In any case, previous analyses of CMB data allowing 
for a varying a || ^, |o| have revealed some interesting 
degeneracies between a and other cosmological param- 
eters, such as £lb or Hq. On the other hand, a recent 
'brute-force' exploration of a particular sector of param- 
eter space (including quintessence models) seems to 
claim results on degeneracies between the various param- 
eters[§|,@. 

While the two approaches are not really comparable 
( [pd| being rather more simplistic, as it uses no actual 
data and has somewhat unclear criteria for the presence 
of a degeneracy), this discrepancy begs the question of 
whether the degeneracies found in [|[ [| [l(J are real 'phys- 
ical' and fundamental degeneracies, which will remain 
at some level, no matter how much more accurate data 
one can get, or if they are simply degeneracies in the 
data, which won't necessarily be there in other (better) 
datasets. And a related question is, of course, assum- 
ing that the degeneracies are significant, how can one get 
around them. We will address these issues in the follow- 
ing sections. 



III. CURRENT CMB AND LARGE-SCALE 
STRUCTURE CONSTRAINTS 

Here we present an up-to-date analysis of the Cosmic 
Microwave Background constraints on varying a as well 
as, for the first time, an analysis of its effects on the 
large-scale structure (LSS) power spectrum. 

Even though this may not be entirely obvious, a vary- 
ing a will have an effect on the matter power spectrum. 
The simplest way to understand this is to interpret the 
variation in a as being due to a variation in the speed of 
light c (which one is always free to do pi|]). 

A variation in a affects the matter power spectrum 
to the extent that it changes the horizon size, hence the 
turnover scale in the matter power spectrum. Allow- 
ing for a variation in a, this is not only a function of 
f2 m , Qb and h but of a as well, through the dependence 
of the recombination epoch on a. Therefore varying al- 
pha will produce a change in the turnover point position 
k rec of the matter power spectrum, hence a shift of the 
curve sideways, and therefore a change on the value of 
erg . For example a decrease in a shifts this turnover scale 
to smaller k, hence allowing for a decrease in erg. 

By plotting the transfer functions (generated with a 



modified § version of the CMBFAST || code which 
includes the effects of a varying a) we find that this effect 
is fairly small. For fi TO = 0.3 and £Ib = 0.05 and keeping 
all other cosmological parameters fixed a variation of a 
by say 10% from its standard value produces variations in 
the transfer function which are at most 5% in restricted 
regions of k (the effect on the value of as is even smaller). 

On the other hand a change in a will modify the height 
of the first peak of the CMB power spectrum through the 
(early) ISW effect. This effect also depends on VL m h 2 . 
This illustrates the interplay between a varying a, fl m 
and the value of as- Further effects of a varying a in 
the CMB are a slight change in the position of the first 
peak due to the aforementioned change in the horizon 
size, plus a variation in the high-£ damping (due to the 
finite thickness of the last-scattering surface) which are 
also dependent on a number of cosmological parameters 
other than a. 

It should be emphasized that although these CMB and 
LSS constraints are in some sense complementary, and 
can help break degeneracies by determining other cosmo- 
logical parameters, they certainly can not be blindly com- 
bined together, since the range of cosmological epochs (or 
redshifts) to which they are sensitive is somewhat differ- 
ent. 



A. CMB data analysis 

We compare the recent CMB observations with a set 
of flat models with parameters sampled as follows (the 
value in brackets is the step size): 

n m = o.i. ..(o.i)... i.o (6) 

Q b = tt m - ttcdm = 0.009 . . . (0.003) . . . 0.036 (7) 

— = 0.80... (0.01)... 1.10 (8) 
a 

h = 0.40... (0.05)... 0.90 (9) 

n s = 0.70 ... (0.05) ... 1.30 . (10) 

We rescale the amplitude of fluctuations by a pre- factor 
do, in units of C? OBE , with 0.50 < C w < 1.40. Fi- 
nally, we assumed a negligible re-ionization and an op- 
tical depth t c ~ 0. This is in agreement with recent 
estimates on the redshift of re-ionization z re ~ 6 ± 1 (see 
e.g. fH). 

The theoretical models are computed using a modi- 
fied version of the publicly available CMBFAST pro- 
gram pj|, accounting for the effects of a varying a, and 
are compared with the recent BOOMERanG-98, DASI 
and MAXIMA- 1 results. The power spectra from these 
experiments were estimated in 19, 9 and 13 bins respec- 
tively, spanning the range 25 < I < 1150. 
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For the DASI and MAXIMA-I experiment we use 
the publicly available correlation matrices and window 
functions. For the BOOMERanG experiment we as- 
sign a constant value for the spectrum in each bin 
£(£+l)Cg/27r = Cb, we approximate the signal Cb inside 
the bin to be a Gaussian variable and we consider ~ 10% 
correlations between contiguous bins. The likelihood for 
a given cosmological model is then defined by 



21n£ = [Cf - C b x )M B b'{C^, - C B X ,) 



th 



(11) 



where Mg b' is the Gaussian curvature of the likelihood 
matrix at the peak. We consider 10%, 4% and 5% Gaus- 
sian distributed calibration errors for the BOOMERanG- 
98, DASI and MAXIMA-1 experiments respectively and 
we included the beam uncertainties by the analytical 
marginalization method presented in |p7[ . We also in- 
clude the COBE data using Lloyd Knox's RADPack 
packages. 



B. LSS data analysis 

In what follows, we will add to the CMB data the real- 
space power spectrum of galaxies in the 2dF 100k galaxy 
redshift survey using the data and window functions of 
the analysis of p8[ . 

To compute Lr^ F we, evaluate pi — P(ki), where P(k) 
is the theoretical matter power spectrum and fc$ are the 
49 k- values of the measurements in ^8|. Therefore we 
have 



2ln£ 



2dF 



E 



[Pi - (WpUVdPi 



(12) 



where Pi and dPi are the measurements and error bars 
in |2f| and W is the reported 27 x 49 window matrix. 
We restricted the analysis to a range of scales where the 
fluctuations are assumed to be in the linear regime (k < 
Q.Q2h" 1 Mpc). When combining with the CMB data, we 
marginalize over a bias b considered as additional free 
parameter. 

We will also include information on a$ , the RMS mass 
fluctuation in spheres of 8h~ 1 Mpc, obtained from local 
cluster number counts. There is presently no consensus 
on the correct value of this observable, mainly because 
of systematics in the calibration between cluster virial 
mass and temperature. For convenience of analysis, we 
consider 2 values: an high value ~ ffi„ 6 crg = 0.50 ± 0.05 
in agreement with the results of |p9] , pOf and a lower one, 
~ n™a 8 = 0.40 ± 0.05 following the analysis of ^ ||. 

We attribute a likelihood to each value of 8a /a by 
marginalizing over the nuisance parameters. We then 
define our 68% (95%), confidence levels to be where the 
integral of the likelihood is 0.16 (0.025) and 0.84 (0.975) 
of the total value (see e.g. pjjfl). 
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FIG. 1: The current combined CMB constraints on the u>b 
a/ao, n s — a/ao and h — a/ao planes. 
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TABLE I: Current cosmological constraints on the variation 
of the fine-structure constant (marginalizing over other pa- 
rameters) for various different priors. 



Prior 


a/ao (95 % c.l.) 


h = 0.65 ±0.2 


o.95lg:g? 


BBN io b = 0.02 ± 0.002 


0.96lH 6 5 


HST h = 0.71 ± 0.08 


0-98±g : gg 


SN-Ia 


Q5+ 07 


0£iV 8 = 0.50 ± 0.05 


1 01+ 004 

1 - ui -0.04 


f^ 6 cr 8 = 0.40 ± 0.05 


o-99±8:SS 


2dF 





C. Results 

In a previous work || we produced likelihood con- 
tours in the fib h 2 — a/ao plane by analyzing the recent 
BOOMERanG and DASI CMB datasets and by includ- 
ing two priors to the analysis: flatness and h = 0.65±0.2. 
Our results were consistent with the baryon abundance 
obtained from big-bang nucleosynthesis (BBN) and we 
constrained variations in a at z ~ 1000 at level of about 
10%. In Fig. |l| we plot constraints on this plane, as well 
as on the n s — a/ao and h — a/ao planes with a similar 
analysis, but including this time the MAXIM A-I dataset. 

From these results we can see that the inclusion of the 
MAXIMA-I data doesn't significantly change our pre- 
vious constraints. In fact, even if the analysis of the 
MAXIMA-I data alone suggests a higher value of the 
baryon fraction (fl^h 2 ~ 0.030±0.005 [M), the combined 
analysis with DASI and BOOMERanG still suggests a 
low value of fl^h 2 ~ 0.023. This result is in agreement 
with previous analysis, e.g. | p8| . 

From last two panels of Fig. [j] we also see that, in 
the set of models we are considering, there is a clear 
correlation between variations in a and changes in the 
scalar spectral index n s and the Hubble parameter h. 
We will see in the next section that this degeneracy can 
be broken by the future and more accurate measurements 
from satellite experiments like MAP J35| or the Planck 
Surveyor |3l| . 

However, since variations in n s and h affect the shape, 
position and amplitude of the matter power spectrum 
(irrespective of changes in a), we can in principle use the 
data from galaxy clustering and local cluster abundances 
in order to break these degeneracies and infer stronger 
constraints on a. Regarding the RMS amplitude of mass 
fluctuations on scales of 8h~ 1 Alpc, the two effects are 
competitive: an increase in ns will increase erg, while 
lowering h decreases it. As we can see from the first 
panel of Fig. ^, the effect from h is stronger and the final 
result is that a decrease in a generally allows for lower 
amplitudes of as- 

Furthermore, as we can see from Fig. |J a certain de- 
gree of degeneracy is present in the CMB data between a 




FIG. 2: The current combined CMB constraints on the erg — 
a/ao and Q. m h — a/ao planes. The latter panel also shows 
the results of combining the CMB and 2dF datasets. 



and the shape parameter fl m h. This degeneracy can be 
optimally broken by incorporating the data from the 2dF 
galaxy survey. As we can see in the bottom panel, includ- 
ing the 2dF data shrinks the contours around a/ao ~ 1 
and Q m h ~ 0.2, In other words, there is a clear distinc- 
tion between the CMB and LSS data: while for the CMB 
data a negative variation of a is preferred, the opposite 
happens for the LSS data. When the two datasets are 
combined, the best-fit model, in the £l m h — a/ao plane, 
is quite close to the standard one. 

We should emphasize at this point that in doing this 
we are not combining direct constraints on the parame- 
ter a itself, obtained through both methods, to obtain a 
tighter constraint. As mentioned above this can not be 
done, since the CMB and LSS analyses are sensitive to 
the values of a at different redshift ranges, so there is no 
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reason why these values should be the same. Addition- 
ally, there is no well-motivated theory that could relate 
such variations at different cosmological epochs. All that 
one could do at this stage would be to assume some toy 
model where a certain behaviour would occur, but this 
would mean introducing various additional parameters, 
thus weakening the analysis. Hence we chose not to pur- 
sue this path and leave the analysis as model-independent 
as possible. 

What we are doing is using additional information 
(which is also sensitive to a) to better constrain other 
parameters in the underlying cosmological model, such 
as n s , h and the densities of various matter components, 
which we can reliably assume are unchanged throughout 
the cosmological epochs in question. In other words, we 
are simply selecting more stringent priors for our analy- 
sis, in a self-consistent way. 

The constraints obtained by combined analysis are re- 
ported in Table Q. Our main result from this table is 
that, as one would expect, when constraints from other 
and independent cosmological datasets are included in 
the CMB analysis, the constraints on variations on a be- 
come significantly stronger. 



IV. FISHER MATRIX ANALYSIS 

The precision with which the forthcoming satellite ex- 
periments MAP H| and Planck Q will be able to de- 
termine variations in a can be readily estimated with a 
Fisher Matrix Analysis (FMA). Some authors have al- 
ready performed such an analysis in the past ]37], |38[ : 
however, their analysis was based on a different set of 
cosmological parameters and assumed cosmic variance 
limited measurements. In our FMA we also take into ac- 
count the expected performance of the MAP and Planck 
satellites and we make use of a cosmological parameters 
set which is well adapted for limiting numerical inaccu- 
racies. Furthermore, the FMA can provide useful insight 
into the degeneracies among different parameters, with 
minimal computational effort. 



A. Analysis Setup 

We characterize the cosmological model by a 7 dimen- 
sional parameter set, given by 



© = [a, u b , uj m , lu Ai TZ, n s , Q), 



(13) 



where u>b = Qbh 2 is the physical baryonic density, tu m = 
(^cdm + Qb)h 2 the energy density in matter and to a the 
energy density due to a cosmological constant. Here h 
denotes the Hubble parameter today, Ho = lOO/i km 
s" 1 Mpc -1 . The quantity TZ = £ rct /£ is the 'shift' pa- 
rameter (see |3£| Q and references therein) , which gives 
the position of the acoustic peaks with respect to a flat, 
fl A = reference model, n s is the scalar spectral index 



and Q =< 1(1 + l)Ct > 1 / 2 denotes the overall normaliza- 
tion, where the mean is taken over the multipole range 
2 < I < 2000. 

The shift parameter TZ depends on f2 m , on the curva- 
ture Q K = 1 — J1a — — f^rad through 
1 



K 



2 1 



VP 



7\ l 



xiv) 



n 



rad 



l 



rad 



(14) 



where Zd cc is the redshift of decoupling, fl rac \ is the energy 
parameter due to radiation (f2 ra d = 4.13 • 10~ 5 /ft. 2 for 
photons and 3 neutrinos) and 



dz 



v = vW 

[r>rad(l + Zf + fi m (l + Zf + a K (l + Zf + 



(15) 
tt A ]- 1/2 . 

The function x{y) depends on the curvature of the uni- 
verse and is y, sin(y) or sinh(j/) for flat, closed or open 
models, respectively. Inclusion of the shift parameter TZ 
into our set of parameters takes into account the geo- 
metrical degeneracy between wa and u) m Kj] ]. With our 
choice of the parameter set, TZ is an independent variable, 
while the Hubble parameter h becomes a dependent one. 

We assume throughout purely adiabatic initial condi- 
tions and we do not allow for a tensor contribution. In the 
FM approach, the likelihood distribution for the param- 
eters is expanded to quadratic order around its maxi- 
mum. We denote this maximum likelihood (ML) point by 
0o and call the corresponding model our "ML model" , 
with parameters lo^ = 0.0200 (hence fib = 0.0473), 
uj m = 0.1267 (hence Q m = 0.3000), uj a = 0.2957 (hence 
fl A = 0.7000 and h = 0.65), TZ = 0.9628^ s = 1.00, 
Q = 1.00. For the value of Zd cc in eq. (^5|) (which is 
weakly dependent on tub and Wtot) we have used the 
fitting formula from jf2|. For the ML model we have 
z dcc = 1115.52. 

TABLE II: Experimental parameters for MAP and Planck 
(nominal mission). Note that the sensitivities are here ex- 
pressed in jiK. 





MAP 


Planck 


v (GHz) 
9 C (arcmin) 
a c (/iK) 

w- 1 ■ 10 15 (K 2 ster) 

£c 


41 61 95 
31.8 21.0 13.8 
19.8 30.0 45.6 
33.6 33.6 33.6 
254 385 586 


100 143 217 
10.7 8.0 5.5 
4.6 5.4 11.7 
0.215 0.158 0.350 
757 1012 1472 


*-max 
/sky 


1000 
0.80 


2000 
0.80 



Proceeding as described in 
Fisher information matrix 



hi, we then calculate the 



1=2 



1 dC t dCj 
AC? 00*50,- 



(16) 
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The quantity ACV is the standard deviation on the esti- 
mate of Ce: 



ACf = 



(2£+l)f sky 



-2\2 



(17) 



The first term is the cosmic variance, arising from the fact 
that we exchange ensemble average with a spatial aver- 
age. The second term takes into account the expected 
error of the experimental apparatus fl3| : 



^ 2 = E" 



(18) 







The sum runs over all channels of the experiment, with 
the inverse weight per solid angle w^ 1 = (a c c )~ 2 and 
l c = v8 In 2/0 c , where a c is the sensitivity (in /iK) and 
6 C is the FWHM of the beam (assuming a Gaussian pro- 
file) for each channel. Furthermore we can neglect the is- 
sues arising from point sources, foreground removal and 
galactic plane contamination assuming that once they 
have been taken into account we are left with a "clean" 
fraction of the sky given by / s ky 

The experimental parameters are summarized in Table 
We use the 3 higher frequency MAP channels and the 
first 3 channels of the Planck High Frequency Instrument 
(HFI) . Adding the 3 higher frequency channels of the HFI 
and the 3 channels of Planck's Low Frequency Instrument 
leaves the expected errors unchanged: therefore they can 
be used for foreground removal, consistency checks, etc, 
leaving the HFI channels for cosmological use. 

For Gaussian fluctuations, the covariance matrix is 
then given by the inverse of the Fisher matrix, C = F^ 1 
p4j . The la error on the parameter 0^ with all other 
parameters marginalised is then given by y/Ca. If all 
other parameters are held fixed to their ML value s, the 
standard deviation on parameter 0j reduces to * / 1 / Fa 
(conditional value). Other cases, in which some of the 
parameters are held fixed and others are being marginal- 
ized over can easily be worked out. 

A case of interest is the one in which all parameters are 
being estimated jointly: then the joint error on parameter 
i is given by the projection on the i-th coordinate axis of 
the 7-dimensional hyper-ellipse which contains a fraction 
7 of the joint likelihood. The equation of the hyper- 
ellipse is 



(0 - O )F(0 - O ) 



qi- 



'7' 



(19) 



where f?i- 7 is the quantile for the probability 1 — 7 for a 
X 2 distribution with 7 degrees of freedom. For 7 = 0.683 



(la c.l.) we have 



.18. 



The FMA assumes that we are expanding the likeli- 
hood function at the right point, ie that the parameters 
values of the true model are in the vicinity of ©o- The 
validity of the results depends on this assumption, as 
well as on the assumption that the a^ m 's are independent 
Gaussian random variables. If the FM predicted errors 
are small enough, the method is self-consistent and we 



can expect the FM prediction to reproduce in a correct 
way the exact behaviour. This is indeed the case for 
the present analysis, with the notable exception of wa, 
which suffers from the geometrical degeneracy (see next 
section). 

Special care must be taken in computing the deriva- 
tives of the power spectrum with respect to the cosmo- 
logical parameters. Numerical errors in the spectra can 
lead to larger derivatives, which would artificially break 
degeneracies among parameters. In the present work we 
implement double-sided derivatives, which diminish the 
truncation error from second order to third order terms. 
The choice of the step size is a trade-off between trun- 
cation error and numerical inaccuracy dominated cases. 
For an estimated numerical precision of the computed 
models of order 10 , the step size should be approx- 
imately 5% of the parameter value p5| ]. It turns out 
that for derivatives in direction a and n s the step size 
can be chosen to be as small as 0.1%. As for the other 
parameters, the accuracy is limited by the fact that dif- 
ferentiating around a flat model requires computing open 
and closed models, which are calculated using different 
numerical techniques. The relative numerical noise is 
therefore much larger. After several tests, we chose step 
sizes varying from 1% to 5% for u>b, u) m , lua, and 1Z. This 
choice gives derivatives with an accuracy of about 0.5%. 
The derivatives with respect to Q are exact, being the 
power spectrum itself. 



B. Analysis results 



1. FMA forecast 



Table III summarizes the results of our FMA. MAP 



will be able to constrain variations in a at the time of 
last scattering to within 2% (la, all others marginalised). 
This corresponds to an improvement of a factor of 3 
relative to the limits presented in the previous section. 
Planck will narrow it down to about half a percent. If 
all other parameters are supposed to be known and fixed 
to their ML value, then a factor of 10 is to be gained in 
the accuracy of a (compare the columns labelled "fixed" 
in table III). However, if all parameters are being esti- 
mated jointly, the accuracy on variations in a will not go 
beyond 1%, even for Planck (column "joint"). 

The parameters lo^^TZ and n s suffer from partial de- 
generacies with a, which are discussed in more detail 
in the next section. This is only partially reflected in 
the marginalized errors of table HI . Correlations among 



the parameters play an important role: within the limit 
of the quadratic order approximation, they are fully de- 
scribed by the FM. 

The geometrical degeneracy limits the accuracy on u> m 
and wa ■ The degeneracy is so severe that the error on co\ 
is very unsensitive to the experimental details. From the 
FMA point of view, this happens because the derivative 
of the spectrum with respect to uj\ vanishes for i > 50. 
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Therefore probing higher multipoles does not help for the 
purpose of better constraining the cosmological constant. 
We emphasize once more that such a large error cannot 
be trusted to be accurate in any respect: it just signals 
a very large inaccuracy in u>\. The errors on all other 
parameters, however, are small enough to justify the self- 
consistency of the FMA approach. 

The power of an experiment can be roughly assessed 
by looking at the eigenvalues Xi and eigenvectors uW of 
its FM: the error along the direction in parameter space 

defined by (principal direction) is proportional to 

1 / 9 

A t x '~. But we are interested in determining the errors on 
the physical parameters rather then on their linear com- 
binations along the principal directions. Therefore in the 
ideal case we want the principal directions to be as much 
aligned as possible to the coordinate system defined by 
the physical parameters. We display in Table IV eigen- 
values and eigenvectors of the FM for MAP and Planck. 
Planck's errors, as measured by the inverse square root 
of the eigenvalues, are smaller by a factor of about 4 on 
average. For 6 of the 7 eigenvectors Planck also obtains 
a better alignment of the principal directions with the 
axis of the physical parameters. This is established by 
comparing the ratios between the largest (marked with 
an asterisk in Table IV) and the second largest (marked 
with a dagger) cosmological parameters' contribution to 
the principal directions. This is of course in a slightly 
different form the statement that Planck will measure 
the cosmological parameters with less correlations among 
them. 



2. Degeneracies with other parameters 

In previous work || some of the present authors ob- 
served a degeneracy in the Boomerang-98 and Maxima- 1 
data between a and u>b- This degeneracy also shows up 
in the present analysis (see Fig. flL top panel). The ques- 
tion we ask is: is this a fundamental degeneracy, or is it 
only in the data? 

We have performed a FMA with experimental param- 
eters chosen as to mimic a Boomerang-type experiment 
(/ sky = 0.02, a = 15 /iK, 6 = 9.2', £ max = 1000). If the 
degeneracy is due to the limited precision of the present- 
day experimental data, we expect the degeneracy to dis- 
appear as we move from Boomerang, to MAP, to Planck. 
Fig. H (top panel) shows la joint confidence curves (all 
other parameters marginalized) in the a/aQ—Ub plane for 
the FM simulated Boomerang, MAP and Planck (from 
the outside to the center, respectively). The curve for 
Boomerang is to be compared with the la contour of the 
data analysis (Figure [l], top panel). Although the FM 
ellipse is centred by construction at the ML model value, 
it is in qualitative agreement with the result of the data 
analysis. As we move to MAP, the degeneracy shrinks 
but is still there: only higher multipole measurements 
from Planck can break it. 

The same behaviour is observed in the a/ctQ — n s plane 
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FIG. 3: Ellipses containing 68% (lcr) of joint confidence (all 
other parameters marginalized) in the a/cto — Uj,, ct/oto — n s 
and a/ao—TZ planes for the FM simulated Boomerang, MAP 
and Planck. 
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TABLE III: Fisher matrix analysis results: expected la errors for the MAP and Planck satellites. The column marg. gives the 
error with all other parameters being marginalized over; in the column fixed the other parameters are held fixed at their ML 
value; in the column joint all parameters are being estimated jointly. 



Quantity 


la errors (%) 






MAP 






Planck HFI 






marg. 


fixed 


joint 


marg. 


fixed 


joint 


a 


2.24 


0.13 


6.39 


0.41 


0.02 


1.16 


LOb 


5.11 


1.12 


14.61 


0.98 


0.31 


2.79 


LO m 


5.26 


1.97 


15.04 


2.30 


0.44 


6.59 


L0 A 


97.81 


89.62 


279.74 


95.17 


89.55 


272.18 


TZ 


3.73 


0.20 


10.67 


0.57 


0.03 


1.64 


n s 


1.79 


0.52 


5.12 


1.19 


0.13 


3.42 


Q 


1.19 


0.36 


3.41 


0.19 


0.10 


0.54 



TABLE IV: In the lines we display the components of the eigenvectors of the FM for MAP and Planck. The quantity l/^/\l 
is proportional to the error along the principal direction u' 1 '. For each principal direction, an asterisk marks the largest 
cosmological parameters contribution, a dagger the second largest. Not only Planck has errors smaller by a factor of about 4 
on average, but also the alignment of the principal directions with the axis denned by the physical parameter is better than 
MAP in 6 cases out of 7. 



MAP 


Direction i 




LOb 




LOA 


n s 


Q 


TZ 


a 


1 


2.23E-04 


9.9540E-01* 


-9.7492E-03 


-1.8938E-05 


-1.6970E-02 


-4.2258E-03 


-6.8163E-02t 


6.4257E-02 


2 


1.12E-03 


-7.8215E-02 


2.8225E-01 


1.6168E-05 


1.6425E-01 


-1.4652E-01 


-4.4305E-01f 


8.1822E-01* 


3 


2.43E-03 


4.0471E-02 


7.2825E-01* 


2.0996E-04 


1.4674E-01 


-5.6711E-01t 


2.6310E-01 


-2.3590E-01 


4 


8.63E-03 


3.9968E-03 


6.1138E-01t 


8.1457E-03 


-3.7196E-01 


6.4228E-01* 


-2.3162E-01 


-1.4702E-01 


5 


4.45E-02 


-1.9844E-02 


-6.8856E-02 


-1.6683E-02 


2.5892E-01 


-1.8786E-01 


-8.0247E-01* 


-4.9829E-01f 


6 


1.38E-02 


3.1792E-02 


1.0650E-01 


-1.5423E-02 


8.6336E-01* 


4.5726E-01t 


1.7932E-01 


-2.8026E-02 


7 


2.89E-01 


2.0351E-04 


-4.6453E-03 


9.9971E-01* 


2.0638E-02t 


-1.1925E-03 


-8.7876E-03 


-7.5125E-03 


Planck 


Direction i 


1/VAi 


LOb 


LOm 


LOA 


n s 


Q 


TZ 


a 


1 


5.81E-05 


9.3678E-01* 


8.2153E-03 


-1.2417E-06 


9.2222E-03 


3.0060E-03 


-1.8611E-01 


2.9604E-01t 


2 


2.87E-04 


-3.2830E-01 


2.5649E-01 


9.3971E-06 


1.5805E-01 


2.0156E-01 


-3.8542E-01t 


7.8248E-01* 


3 


5.46E-04 


1.1323E-01 


8.4059E-01* 


4.5325E-06 


2.7521E-01 


2.5135E-01 


3.2628E-01t 


-1.8765E-01 


4 


1.94E-03 


3.8999E-02 


-3.1741E-01t 


-9.3635E-05 


-3.8346E-02 


9.4091E-01* 


6.4234E-02 


-8.2576E-02 


5 


2.86E-03 


-1.5404E-02 


2.9918E-01 


7.3411E-04 


-3.9324E-01 


9.9942E-02 


-7.5339E-01* 


-4.2194E-011 


6 


1.30E-02 


8.7278E-03 


-1.9309E-01 


-1.6271E-02 


8.6190E-01* 


-2.9783E-02 


-3.7230E-01| 


-2.8285E-01 


7 


2.81E-01 


1.6073E-04 


-3.3978E-03 


9.9987E-01* 


1.4308E-02t 


-4.7295E-04 


-5.4974E-03 


-4.3069E-03 



(Fig. H| middle panel; compare with Fig. [[], middle 
panel). Again, the observed degeneracy between a and 
n s is clearly revealed by the FMA for Boomerang. 

In the bottom panel of Fig. [| we investigate the im- 
portant degeneracy between TZ and a/a^. These two pa- 
rameters are very highly correlated (correlation ss 0.99 
for all experiments), because an increase of a displaces 
the acoustic peaks to higher multipoles. This effect is 
mainly due to the increased redshift of last scattering 
H Q . On the contrary, an increase of TZ shifts the peaks 
toward smaller i values, because of the change in the an- 
gular diameter distance relation. However, an increase in 
a also produces a decrease in the damping at high mul- 



tipoles, which can be used to break the degeneracy, as it 
is the case for Planck. This is reflected in the different 
amplitude for the two derivatives, which are otherwise 
perfectly in phase, as can be seen in Fig. |j. This degen- 
eracy can clearly be identified because of our choice of 
the parameter set ©, which includes the shift parameter 
rather than the Hubble or curvature parameters. This 
emphasizes the importance of a correct choice of the pa- 
rameter set in the context of a FMA. 

Comparison to previous works is only partially possi- 
ble, because of the differences in the analysis discussed 
above. Our detailed analysis confirms however the con- 
clusions in refs. [B7l p8|, which found that a cosmic vari- 
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FIG. 4: Derivatives with respect to a and 1Z. We plot 
—dCi/da to facilitate the comparison with dCe/dTZ. The 
two derivatives are perfectly in phase: this is responsible for 
the degeneracy between the corresponding parameters. Only 
the different amplitudes allow an experiment which maps suf- 
ficiently high multipoles with high accuracy to distinguish 
between them. 



ance limited experiment could obtain a precision on a 
of order 10 -3 — 10~ 2 . We have also shown that there 
is much to be gained from using prior knowledge about 
the other parameters in the determination of a via CMB 
measurements. The improvement in accuracy is about a 
factor 50 for both MAP and Planck. 



V. CONCLUSIONS 

We have provided an up-to-date analysis of the effects 
of a varying fine-structure constant a in the CMB, focus- 
ing on the issue of the degeneracies with other cosmolog- 
ical parameters, and of how these can be broken. 

We have shown that the currently available data is con- 
sistent with no variation of a from the epoch of recom- 
bination to the present day, though interestingly enough 
the CMB and LSS datasets seem to prefer, on their own, 
variations of a with opposite signs. Whether or not this 
statement has any physical relevance (beyond the results 
of the statistical analysis) is something that remains to 
be investigated in more detail. In any case, any such 
(relative) variation is constrained to be less than about 
4%, so a best-fit or 'concordance' model with a exactly 
constant will require, at most, some slight deviations of 
other cosmological parameters from the 'standard' values 
obtained from analyses which don't allow for a variation. 

On the other hand, the prospects for the future are 
definitely bright. In the short term, the imminent VSA 
and CBI data should be able to provide some improve- 



ment on the current results. The dataset we have 
used (BOOMERanG, MAXIMA and DASI) all have the 
common feature that their error bars are smallest for 
data points around the first Doppler peak and larger 
for smallest angular scales. Now, as we have explicitly 
shown above (and was already suggested in |9) ) , the first 
Doppler peak is not a vary accurate 'a-meter', due to 
the degeneracy with the shift parameter. Thus datasets 
in which points around the first peak will somewhat dom- 
inate the statistical analysis are not optimal for a esti- 
mation. In this regard, VSA and CBI should be useful 
because they can provide a significant number of data 
points on small angular scales with relatively small error 
bars, hence minimizing this problem. 

In the longer term, the forthcoming satellite experi- 
ments will provide a dramatic improvement on these re- 
sults. We have performed a Fisher Matrix Analysis using 
a well adapted parameter set and realistic experimental 
characteristics for the upcoming MAP and Planck satel- 
lite missions. The results of our forecast are that MAP 
and Planck will be able to constrain variations in a within 
2.2% and 0.4%, respectively (la c.L, all other parame- 
ters marginalized). If all parameters are being estimated 
simultaneously, then this limits increase to about 6.4% 
and 1.1%, respectively. The analysis of the presently ob- 
served degeneracies between a and n s comes to the 
conclusion that measurement of higher multipoles will 
allow to break it. We have also identified an important 
degeneracy between a and the shift parameter. 

To conclude, we have provided a thorough analysis 
of the effects of cosmological parameter degeneracies in 
CMB measurements of the fine-structure constant a, and 
quantified the importance of these degeneracies. We have 
also explicitly discussed two ways in which these degen- 
eracies can be circumvented, namely acquiring better 
data (the easy solution, at least from the theorists' point 
of view) or combining the CMB data with other cosmo- 
logical datasets which can provide constraints on other 
cosmological parameters (the 'brute-force' solution) 

In a follow-up paper, we will discuss a third way in 
which these degeneracies can be lifted, namely including 
CMB polarization data (l2) (the more elegant solution 
in principle, though it's yet to be realized in practice). 
These tools, together with other measurements coming 
from BBN || and quasar and related data ]l], || offer 
the exciting prospect of being able to map the value 
of a at very many different cosmological epochs, which 
would allow us to impose very tight constraints on higher- 
dimensional models where these variations are ubiqui- 
tous. 
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